Strength of reduced two-body spin-orbit interaction from chiral three-nucleon force 



M. Kohno 

Physics Division, Kyushu Dental College, Kitakyushu 803-8580, Japan 

The contribution of a chiral three-nucleon force to the strength of an effective spin-orbit coupling 
is estimated. We first construct a reduced two-body interaction by folding one-nucleon degrees of 
freedom of the three-nucleon force in nuclear matter. The spin-orbit strength is evaluated by a 
Scheerbaum factor obtained by the G-matrix calculation in nuclear matter with the two-nucleon 
interaction plus the reduced two-nucleon interaction. The problem of the insufficiency of modern 
realistic two-nucleon interactions to account for the empirical spin-orbit strength is resolved. It is 
also indicated that the spin-orbit coupling is weaker in the neutron-rich environment. Because the 
spin-orbit component from the three-nucleon force is determined by the low-energy constants fixed 
in the two-nucleon sector, there is little uncertainty in the present estimation. 
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Spin-orbit field in atomic nuclei is essential to repro- 
duce well-established single-particle shell structure. The 
empirical strength of the spin-orbit potential, however, 
has not been fully understood on the basis of the realis- 
tic nucleon-nucleon force. The possible role of intermedi- 
ate isobar A-excitation to the nuclear spin-orbit field was 
considered in parallel with the construction of the two- 
pion-exchange three-nucleon force (3NF) by Fujita and 
Miyazawa [lj. The problem was reinvestigated in the 
early 1980s [2, S] to search for the additional spin-orbit 
strength. Later, the Illinois group showed Q| that their 
3NF makes a substantial contribution to the spin-orbit 
splitting in 15 N. 

Kaiser and his collaborators investigated, in their sev- 
eral papers (EH?}, the nuclear spin-orbit coupling in the 
framework of chiral perturbation theory. The large con- 
tributions generated by iterated one-pion exchange and 
the 3NF almost cancel each other 0, 0] , and the short- 
range spin-orbit strength in the form of the effective four- 
nucleon contact-coupling deduced from realistic nucleon- 
nucleon interactions accounts well Q for the empirical 
one. Because the contact-interaction in the chiral per- 
turbation, however, is still needed to be regulated for the 
application to low-energy nuclear structure calculations, 
and their arguments for various contributions seem not 
to be fully unified, it is worthwhile to analyze the effec- 
tive strength of the spin-orbit coupling by applying the 
established microscopic theory, namely the lowest-order 
Brueckner theory, to the two-nucleon and three-nucleon 
interactions in the chiral effective field theory (Ch-EFT). 

The Thomas form of an average single-particle spin- 
orbit potential has been used to describe nucleon spin- 
orbit coupling: 

Ui'-^i-cr, (1) 
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where the radial function p(r) is a nucleon total den- 
sity distribution. The relation of the strength U® s to 
a two-body effective spin-orbit interaction was derived 
by Scheerbaum Q. By defining the constant Bs(q) for 
the triplet odd component of the effective two-body spin- 



orbit interaction vf® (r) 
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with ji being a spherical Besscl function, the single- 
particle spin-orbit potential for spin-saturated nuclei may 
be written as 
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where r specifics cither a proton or neutron. We re- 
fer to Bs(q) Scheerbaum factor, which is differ- 
ent from the original constant in Ref. @ by a factor 
of — Scheerbaum prescribed q ~ 0.7 fm _1 on the 
basis of the wavelength of the density distribution. We 
employ this prescription. If we assume a naive relation 
Pp( r ) = Pn{r) = hp(r), we recover the Thomas form, Eq. 
((T|), with U° s = jBs{q). It has also been customary to 
use a 5-type two-body spin-orbit interaction 

iW(o x + <r 2 ) • (V r x 5{r)V r ) (4) 

in nuclear Hartree-Fock calculations using 5-type Skyrme 
interactions 0, [13] and even with finite range effective 
forces, e.g., the Gogny force [ll[. This two-body force 
provides a single-particle spin-orbit potential: 
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Therefore, the strength W may be identified as the 
Scheerbaum factor Bs(q). The empirical value of W 
is around 120 MeV-fm 5 in various nuclear Hartree- 
Fock calculations. As will be shown below, the mod- 
ern nucleon-nucleon interactions underestimate the spin- 
orbit strength by about 25 %. 

Applying Schcerbaum's formulation to the 
momentum-space G-matrix calculation in nuclear 
matter with the Fermi momentum fop, we obtain the 
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corresponding spin-orbit strength as follows |12j : 
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xW(q, q){(J + 2)G J J +hlJ+1 (q) + G J J AJ (q) 

-(J-l)G3fJ_ MJ _i(9)}. (6) 

Here, q m ax = \{kp + q) and the weight factor W(q,q) is 
,_ , f 9{k F -q) for 0<g< 
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where #(fc_F — g) is a step function. In Eq. (6), G'[[, 
is the abbreviation of the momentum-space diagonal G- 
matrix element in the spin-triplet channel with the total 
isospin T, total spin J, and orbital momenta I' and I. 

Calculating Bs(q) in the lowest-order Brueckner the- 
ory with the continuous prescription for intermediate 



spectra, as presented below explicitly in Table 1, modern 
two-body nucleon-nucleon potentials are found to give 
smaller values of around 90 Mev-fm 5 compared with the 
empirical one. As has been well known that LOBT cal- 
culations in symmetric nuclear matter with realistic two- 
nucleon force do not reproduce correct saturation prop- 
erty. However, in most case, calculated energies at the 
empirical saturation point kp = 1.35 fm _1 are close to 
the empirical energy of about —16 MeV. This suggests 
that G matrices provide basic information on the effective 
nucleon-nucleon interaction in the nuclear medium, by 
incorporating important short-range correlations, Pauli 
effects and dispersion effects. 



Now we consider the contribution of the 3NF. In this 
article, we estimate it in a two-step procedure. First, 
the 3NF «i23 defined in momentum space is reduced to 
an effective two-nucleon interaction 1*12(3) by folding one- 
nucleon degrees of freedom: 
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Here, we have to assume that remaining two nuclcons are 
in the center-of-mass frame, namely k' x + k' 2 = fei + &2 = 
0. The density-dependent effective two-nucleon interac- 
tion as the effect of the 3NF has been commonly in- 
troduced in the literature fl~3l - [l5l |. Note that the suf- 
fix A means an antisymmetrized matrix element; namely 
\ab}A = \ab—ba) and \abc) a = \abc—acb+bca—bac+cab— 
cba) , and the factor | in Eq. © is an additional statisti- 
cal one. This statistical factor has been often slipped in 
the literature. The recent derivation of the effective two- 
body interaction from the Ch-EFT 3NF by Holt, Kaiser 
and Weise [l5| also seems not to be an exception. If an 
adjustable strength is introduced, the statistical factor 
may be hidden in the fitting procedure. In our case of 
using the Ch-EFT 3NF, the low-energy constants except 
for cd and ce are fixed. Although there may be a room 
to adjust Cd and ce, the contributions to the energy from 
these terms are rather small, if they are in a reasonable 
range. In addition, cd and ce do not contribute to the 
reduce two-nucleon spin-orbit interaction. By comparing 
the nuclear matter energy directly calculated from V123 
and that by the reduced 1*12(3), the error due to this ap- 
proximation can be checked to be less than 10 %, if we 
calculate Born energy without including a form factor. 



To explain the procedure of obtaining 1*12(3) more ex- 
plicitly, we write the reduced spin-orbit component orig- 



inating from the c\ term of the Ch-EFT 3NF: 
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where gA = 1.29, / ff = 92.4 MeV, is a pion mass, 
and q { is a momentum transfer of the i-th nucleon. The 
momentum transfer of the third nucleon k is dictated by 
the relation q k = —q i — q^. The folding of the 3NF by 
one nucleon is carried out without incorporating a three- 
body form factor. A form factor is later introduced on 
the two-body level. The folding in symmetric nuclear 
matter with the Fermi momentum kp gives, besides the 
central and tensor components, the following spin-orbit 
term: 
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When carrying out the folding in pure neutron matter, 
the restriction of the isotopic spin brings about an addi- 
tional factor of -|. 

The partial- wave decomposition of the above spin-orbit 
term becomes 
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for the orbital and total angular momenta t and J. The 
functions Q l w and Wf s are defined by 

QwMM) = f F dk 3 Qt(x')Q e (x), (12) 

x {KQt{x)(Qt-i(m') - Qi+i(x')) 
+h 1 Q l {x'){Q l - 1 {x) - Qi+i(x))} , (13) 



where Qi{x) is a Legendre function of the second kind, 
and x' = - a 4j%3 fc^"" and x = fc3 ^^™ T , respectively. 



The spin-orbit component arises also from the c 3 term 
of the Ch-EFT 3NF. This case, in addition to the re- 
placement of the coupling constant, an additional factor 
(ki — k 3 ) ■ (k 3 — fei) appears in the denominator in Eq. 
(|10[) . The partial- wave decomposition reads 
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where the new functions F (k) and Fi(fc) are defined by 
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Adding the reduced two-nuclcon interaction to the Ch- 
EFT two-nucleon interaction, we repeat the LOBT G- 
matrix calculation. Although explicit expressions are not 
shown in this Letter except for the spin-orbit part, we 
include all central, tensor and spin-orbit components of 
the reduced interaction 1*12(3)- The form factor in a func- 
tional form of f{k[M) = exp{-[(fci/A) 4 + (fci/A) 4 ]} is 
introduced for 1*12(3) with the cut-off mass A =550 MeV. 
We use the low-energy constants fixed for the Julich Ch- 
EFT potential by Hebeler et al. [H; c D = -4.381, and 
ce = —1.126. Other constants are c\ = —0.81 GcV -1 , 
C3 = —3.4 GcV" 1 , and C4 = 3.4 GcV -1 . Because the re- 
duction of the 3NF to the two-nucleon force was carried 
out in nuclear matter, 1*12(3) ma y n °t t> e directly applied 
to very light nuclei, such as 3 H and 4 He. 

First, we comment on calculated saturation curves, 
which are given in Fig. 1. Without the contribution 
of the 3NF, the saturation curve attains its minimum 
at larger kp as a function of the Fermi momentum kp 
than the empirical saturation momentum, as has been 
known. Nucleon-nucleon interactions, AV18 [17j, NSC97 
ll, and Julich N 3 LO with the cutoff mass of 550 MeV 
19j give similar saturation curves, and the CD-Bonn po- 
tential poj predicts somewhat deeper binding. For the 
reference of what saturation curve is preferable for nu- 
clear mean filed calculations, we also show the result with 
the Gogny D1S interaction [ill ]. 



k F = 1.35 fm" 1 


AV18 NSC97 CD-B N 3 LO N 3 LO+3NF 


B S (T = 0) 
BsiT = 1) 


2.09 1.9 3.1 2.5 7.0 
86.4 86.7 90.2 84.6 116.2 


Bs{q) 


88.4 88.6 93.3 87.1 123.2 


k F = 1.07 fm" 1 


AV18 NSC97 CD-B N :i LO N 3 LO+3NF 


B S (T = 0) 
BsiT = 1) 


1.4 1.3 2.3 1.6 4.1 
88.1 88.7 92.2 86.5 106.7 


Bs(q) 


89.5 90.0 94.5 88.1 110.8 



TABLE I: B s {q) in the unit of MeV-fm 5 given by Eq. (6) with 
q = 0.7 fm -1 for modern nucleon-nucleon interaction: AV18 
0], NSC97 QJ], CD-Bonn HJ, and Julich N 3 LO Ql]. The 
last entry is the result with including the reduced two-body 
interaction from the Ch-EFT 3NF. 



The thin dotted curve shows the result in which the 
plane wave expectation value of the 3NF 1*123 is added 
to the result of the two-nucleon N 3 LO. The thick dotted 
curve alongside the thin dotted curve is the result with 
the plane wave expectation value of the reduced two- 
nuclcon interaction 1*12(3) • The difference between the 
thin and thick curves is due to the difference of the form- 
factors and the necessary approximation fc 1 +fe 2 = fci + 
k 2 = in Eq. (8). 

The solid curve is the result of the G-matrix calcu- 
lation with including the reduced two-nucleon interac- 
tion, t>i2(3) • Although the energy is seen to be underesti- 
mated by a few MeV, the saturation property is largely 
improved by the repulsive contribution from the three- 
nucleon force. It is not necessary at present to expect 
a perfect agreement with the empirical properties in the 
LOBT calculation in nuclear matter. 
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FIG. 1: Saturation curves in symmetric nuclear matter. 



Now we examine the spin-orbit strength. We tabulate 
values for Bs(q) of Eq. (6) at q = 0.7 fm -1 calculated 
in the LOBT with modern nucleon-nucleon interactions: 
AV18 [nl, NSC97 [H|, CD-Bonn 0, and Jiilich N 3 LO 
(l9l | . The Scheerbaum factors obtained by realistic two- 
nucleon forces are seen to be similar but insufficient to 
explain the strength needed in nuclear mean field calcu- 
lations. Namely only about three-fourths of the empir- 
ically needed strength is accounted for. The two-body 
part of the Ch-EFT, N 3 LO, shows little difference with 
other realistic two-nucleon force. It is also noticed that 
values at kp = 1.07 fm -1 , namely at the half of the 
normal density, change little from those at the normal 
density with kp = 1.35 fin. . It turns out, as the last 
column of Table I shows, that the addition of the reduced 
two-body interaction from the Ch-EFT 3NF bring about 
a good effect to fill the gap, though the 3NF contribu- 
tion is smaller at kp = 1.07 fm -1 . This is in accord 
with the important role of the 3NF to the spin-orbit 
splitting demonstrated in quantum Monte Carlo calcu- 
lations of low-energy neutron-alpha scattering [2lJ. Al- 
though there are ambiguities from the form factor and 
uncertainties inherent in the folding procedure without 
taking into account nucleon-nucleon correlations, no ad- 
ditional adjustable parameter exists, because low-energy 
constants c\ and C3 which contribute solely to the spin- 
orbit strength are determined on the two-nucleon sector. 

As noted after Eq. ([TO]) . the reduced two-body spin- 



orbit term in neutron matter is one-third of that in sym- 
metric nuclear matter. Actual G-matrix calculations us- 
ing the Ch-EFT N 3 LO plus ^12(3) in pure neutron matter 
with k F = 1.35 fm -1 tell that B s (q) values at q = 0.7 
fin -1 are 84.7 and 93.5 MeV-fm 5 without and with the 
reduced two-nucleon interaction ^12(3), respectively. If 
k' F = 1.07 fin -1 is assumed, the corresponding values 
are 87.0 and 94.6 MeV-fm 5 , respectively. Again, the /in- 
dependence is weak. While the spin-orbit strength from 
the two-nucleon force is scarcely different from that in 
symmetric nuclear matter, the additional contribution 
from the three-nucleon force is in fact almost one-third of 
that in symmetric nuclear matter. Thus, the spin-orbit 
strength is expected to be smaller in the neutron-rich en- 
vironment. This seems to be consistent with the trend 
observed in the shell structure near the neutron drip line 
[22| | that a decreasing spin-orbit interaction is preferable 
with increasing neutron excess. 

In summary, we have estimated quantitatively the con- 
tribution of the three-nucleon force of the chiral effective 
field theory to the single-particle spin-orbit strength, us- 
ing the formulation by Scheerbaum Q. We first intro- 
duced the reduced two-body interaction by folding one- 
nucleon degrees of freedom of the 3NF in nuclear matter. 
Making partial- wave expansion of the resulting two-body 
interaction and adding it to the genuine two-nucleon in- 
teraction with including the necessary statistical factor 
of i, wc carried out LOBT G-matrix calculations in in- 
finite matter and evaluated the Scheerbaum factor cor- 
responding to the spin-orbit strength. Because the spin- 
orbit field in the atomic nuclei is fundamentally impor- 
tant as the nuclear magic numbers exhibit, it is impor- 
tant to learn that the inclusion of the 3NF in the chi- 
ral effective field theory can account for the spin-orbit 
strength empirically required for nuclear mean filed cal- 
culations. Because the relevant low-energy constants C\ 
and C3 are determined in the two-nucleon interaction sec- 
tor, there should be little uncertainty for the additional 
spin-orbit strength except for the treatment of the two- 
body form factor. We have also noted that the additional 
spin-orbit strength from the 3NF should be weaker in 
neutron-excess nuclei. 
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